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Abstract: We consider the multi-collinear limit of multi-gluon QCD amplitudes at 
tree level. We use the MHV rules for constructing colour ordered tree amplitudes and 
the general collinear factorization formula to derive timelike splitting functions that 
are valid for specific numbers of negative helicity gluons and an arbitrary number 
of positive helicity gluons (or vice versa). As an example we present new results 
describing the collinear limits of up to six gluons. 
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1. Introduction 


The interpretation of J\f = 4 supersymmetric Yang-Mills theory and QCD as a 
topological string propagating in twistor space [1], has inspired a new and powerful 
framework for computing tree-level and one-loop scattering amplitudes in Yang-Mills 
gauge theory. Notably, two distinct formalisms have been developed for calculations 
of scattering amplitudes in gauge theory - the ‘MHV rules’ of Cachazo, Svrcek and 
Witten (CSW) [2], and the ‘BCF recursion relations’ of Britto, Cachazo, Feng and 
Witten [3, 4], 

In this paper, we wish to exploit these formalisms to examine the singularity 
structure of tree-level amplitudes when many gluons are simultaneously collinear. 
Understanding the infrared singular behaviour of multi-parton amplitudes is a pre¬ 
requisite for computing infrared-finite cross sections at fixed order in perturbation 
theory. In general, when one or more final state particles are either soft or collinear, 
the amplitudes factorise. The first factor in this product is an amplitude depending 
on the remaining hard partons in the process (including any hard partons constructed 
from an ensemble of unresolved partons). The second factor contains all of the sin¬ 
gularities due to the unresolved particles. One of the best known examples of this 
type of factorisation is the limit of tree amplitudes when two particles are collinear. 
This factorisation is universal and can be generalised to any number of loops [5]. 

Both frameworks, the MHV rules and the BCF recursion relations, are remark¬ 
ably powerful in deriving analytic expressions for massless multi-particle tree-level 
amplitudes. At the same time, for the specific purpose of deriving general multi- 
collinear limits, we find the MHV rules approach to be particularly convenient. 

A useful feature of the MHV rules is that it is not required to set reference spinors 
Tj 0 and rja to specific values dictated by kinematics or other reasons. In this way, 
on-shell (gauge-invariant) amplitudes are derived for arbitrary r/’s, i.e. without fixing 
the gauge. By starting from the appropriate colour ordered amplitude and taking 
the collinear limit, the full amplitude factorises into an MHV vertex multiplied by a 
multi-collinear splitting function that depends on the helieities of the collinear gluons. 
Because the MHV vertex is a single factor, the collinear splitting functions have a 
similar structure to MHV amplitudes. Furthermore, the gauge or //-dependence of 
the splitting function drops out. 

One of the main points of our approach is that, in order to derive all required 
splitting functions we do not need to know the full amplitude. Out of the full set 
of MHV-diagrams contributing to the full amplitude, only a subset will contribute 
to the multi-collinear limit. This subset includes only those MHV-diagrams which 
contain an internal propagator which goes on-shell in the multi-collinear limit. In 
other words, the IR singularities in the MHV approach arise entirely from internal 
propagators going on-shell. This observation is specific to the MHV rules method 
and does not apply to the BCF recursive approach. We will see in Section 4.2.3 



that in the BCF picture collinear splitting functions generically receive contributions 
from the full set of allowed BCF diagrams 1 . In view of this, we will employ the MHV 
rules of [2] for setting up the formalism and for derivations of general multi-collinear 
amplitudes. At the same time, various specific examples of multi-collinear splitting 
amplitudes derived in this paper will also be checked in Section 4.2 using the BCF 
recursion relations [3]. 

The basic building blocks of the MHV rules approach [2] are the colour-ordered 
n-point vertices which are connected by scalar propagators. These MHV vertices are 
off-shell continuations of the maximally helicity-violating (MHV) n-gluon scattering 
amplitudes of Parke and Taylor [6, 7]. They contain precisely two negative hclicity 
gluons. Written in terms of spinor inner products [8], they are composed entirely of 
the holomorphic products (i j ) of the right-handed (undotted) spinors, rather than 
their anti-holomorphic partners [ij\, 


A n (l + ,...,p , • • •, n + ) 


_ (pq) 4 _ 

(12) (2 3) • • • (n — 1, n) (n 1) ’ 


( 1 . 1 ) 


where we introduce the common notation {piPj) = (ij) and \piPj\ = [ij]- By 
connecting MHV vertices, amplitudes involving more negative helicity gluons can be 
built up. 

The MHV rules for gluons [2] have been extended to amplitudes with fermions [9]. 
New compact results for tree-level gauge-theory results for non-MHV amplitudes 
involving arbitrary numbers of gluons [10-12], and fermions [9, 13-15] have been 
derived. They have been applied to processes involving external Higgs bosons [16, 17] 
and electroweak bosons [18]. MHV rules for tree amplitudes have further been recast 
in the form of recursive relations [12, 17, 18] which facilitate calculations of higher 
order non-MHV amplitudes in terms of the known lower-order results. In many cases 
new classes of tree amplitudes were derived, and in all cases, numerical agreement 
with previously known amplitudes has been found. 

MHV rules have also been shown to work at the loop-level for supersymmetric 
theories. Building on the earlier work of Bern et al [19, 20], there has been enormous 
progress in computing cut-constructible multi-leg loop amplitudes in J\f = 4 [21-30] 
and A/" = 1 [31-34] supersymmetric gauge theories. Encouraging progress has also 
been made for non-supersymmetric loop amplitudes [35-37]. 

Remarkably, the expressions obtained for the infrared singular parts of A/" = 4 
one-loop amplitudes (which are known to be proportional to tree-level results) were 
found to produce even more compact expressions for gluonic tree amplitudes [30, 38]. 
This observation led to the BCF recursion relations [3, 4] discussed earlier as well 
as extremely compact six-parton amplitudes [39, 40] and expressions for MHV and 
NMHV graviton amplitudes [41, 42], 

lr This is because the required IR poles in the BCF approach arise not only from propagators 
going on-shell, but also from the constituent BCF vertices. 
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The factorisation properties of amplitudes in the infrared play several roles in 
developing higher order perturbative predictions for observable quantities. First, 
a detailed knowledge of the structure of unresolved emission enables phase space 
integrations to be organised such that the infrared singularities due to soft or collinear 
emission can be analytically extracted [45-47]. Second, they enable large logarithmic- 
corrections to be identified and resummed. Third, the collinear limit plays a crucial 
role in the unitarity-based method for loop calculations [19, 20, 48, 49]. 

In general, to compute a cross section at N"LO, one requires detailed knowledge 
of the infrared factorisation functions describing the unresolved configurations for 
n-particles at tree-level, (n — l)-particles at one-loop etc. The universal behaviour in 
the double collinear limit is well known at tree-level (see for example Refs. [50, 51]), 
one-loop [19, 52-56] and at two-loops [57, 58]. Similarly, the triple collinear limit has 
been studied at tree-level [59-62] and, in the case of distinct quarks, at one-loop [63]. 
Finally, the tree-level quadruple gluon collinear limit is derived in Ref. [64], 

Our paper is organised as follows. In Section 2, we briefly review the spinor hclic- 
ity and colour ordered formalism that underpins the MHV rules. Section 3 describes 
the procedure for taking the collinear limit and deriving the splitting functions. We 
write down a general collinear factorization formula, which is valid for specific num¬ 
bers of negative hclicity gluons and an arbitrary number of positive hclicity gluons 
and demonstrate that the gauge dependence explicitly cancels. We find it useful 
to classify our results according to the difference between the number of negative 
hclicity gluons before taking the collinear limit, and the number after. We call this 
difference AM. We provide formulae describing an arbitrary number of gluons for 
AM < 2 in Section 4.1. Specific explicit results for the collinear limits of up to six 
gluons are given in Sec. 4.2.6. We have numerically checked that our results agree 
with the results available in the literature for three and four collinear gluons [64], 
Our findings are summarized in Sec. 5. 

2. Colour ordered amplitudes in the spinor helicity formalism 

Tree-level multi-gluon amplitudes can be decomposed into colour-ordered partial 
amplitudes as 

A n ({ Pi , Ai, a,}) = ig n ~ 2 ^ Tr(TMn ... TM»>) A n (a( l\ ..., n A ")). (2.1) 

cr(E.Sn / 

Here S n /Z n is the group of non-cyclic permutations on n symbols, and j Xj labels 
the momentum pj and helicity A j of the j th gluon, which carries the adjoint repre¬ 
sentation index a*. The T ai are fundamental representation SU(1V C ) colour matrices, 
normalized so that Tr (T a T b ) = 5 ab . The strong coupling constant is a s = g 2 /(47r). 
The MHV rules method of Ref. [2] is used to evaluate only the purely kinematic 
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amplitudes A n . Full amplitudes are then determined uniquely from the kinematic 
part A n , and the known expressions for the colour traces. 

In the spinor helicity formalism [6-8] an on-shell momentum of a massless par¬ 
ticle, = 0, is represented as 

Paa = Pl-i&aa ^a,Xa .) ( 2 - 2 ) 

where A a and Aa are two commuting spinors of positive and negative chirality. Spinor 
inner products are defined by 2 

(A, A') = e ab X a X! b , [A, A'] = -e. h x d \' h , (2.3) 

and a scalar product of two null vectors, p a d = A a Aa and q a a = X' a X' d , becomes 

JW"= -i(A,A'}[A,A'] . (2.4) 

The MHV rules of Ref. [2] were developed for calculating purely gluonic ampli¬ 
tudes at tree level. In this approach all non-MHV n-gluon amplitudes (including 
MHV) are expressed as sums of tree diagrams in an effective scalar perturbation 
theory. The vertices in this theory are the MHV amplitudes of Eq. (1.1) continued 
off-shell as described below, and connected by scalar propagators 1/g 2 . 

When one leg of an MHV vertex is connected by a propagator to a leg of another 
MHV vertex, both legs become internal to the diagram and have to be continued off- 
shell. Off-shell continuation is defined as follows [2]: we pick an arbitrary reference 
spinor r] a and define X a for any internal line carrying momentum g a d by 

Aa = qaaV a • (2.5) 

External lines in a diagram remain on-shell, and for them A is defined in the usual 
way. For the off-shell lines, the same reference spinor rj is used in all diagrams 
contributing to a given amplitude. 

3. The multiple collinear limit 

To End the splitting functions we work with the colour stripped amplitudes A n . For 
these colour ordered amplitudes, it is known that when the collinear particles are not 
adjacent there is no collinear divergence [64], Therefore, without loss of generality, 
we can take particles 1... n collinear. 

The multiple collinear limit is approached when the momenta pi,... ,p n become 
parallel. This implies that all the particle subenergies s %3 = (p t + p 3 ) 2 , with i,j = 

2 Our conventions for spinor helicities follow [1, 2], except that [ij] = — [ij]csw as in ref. [65]. 
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1 ,n, are simultaneously small. We thus introduce a pair of light-like momenta 
P v and i v (P 2 = 0,£ 2 = 0), and we write 


{pi + ■ ■ ■ + PnT 


P v + 


Sl,n r 

H-P' 


Sij = (Pi + ---+ Pjf , 


(3.1) 


where Si jn is the total invariant mass of the system of collincar partons. In the 
collinear limit, the vector P u denotes the collincar direction, and the individual 
collinear momenta are p\ —> z t P u . Here the longitudinal-momentum fractions Zi are 
given by 

£ ' T* <1 

Zi = (3.2) 


and fulhl the constraint Y^i=i z % — 1- To be definite, in the rest of the paper we work 
in the time-like region so that (% >0, 1 > Zi > 0). 


(n + 1) 



(n + 1) 



n 



Figure 1: Factorisation of an IV-point colour ordered amplitude with gluons pi,---,p n 
collinear into splitting function for P —> l,...,n multiplied by an (N — n + l)-point 
amplitude. 


As illustrated in Fig. 1, in the multi-collinear limit an AT-gluon colour ordered 
tree amplitude factorises and can be written as 

A n ( 1 a ‘, .... N J ") split(l Al ,... ,n J " -* P x ) x A N . H+1 ((n + l) A "+‘,..., N X ”,P X ). 

(3.3) 

This labelling of the splitting amplitude split (1- Al ,.—, n Xn —> P A ) differs from the 
usual definition because we use the momentum and hclicity that participates in the 
resultant amplitude P x rather than — P -A . With this choice, it is easier to see how 
the hclicity is conserved in the splitting, i.e. hclicity A 1 ,..., A" is replaced by A. 

There are two different types of collinear limit [2], those that conserve the number 
of negative hclicity gluons between the initial state and the final collinear state, and 
those that do not. 

Only the limits of the type split(l + ,..., n + —> P + ) and split(1 — , 2 + ,..., n + —> 
P~) can contribute to the negative helicity conserving case, and these collinear split¬ 
ting functions are straightforward to derive directly from the simple MHV vertex. 
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All other limits belong to the second class which do not conserve the number 
of negative hclicity gluons, and therefore we classify our results according to the 
difference between the number of negative helicity gluons before taking the collinear 
limit, and the number after, AM. We find that AM corresponds to the order of 
MHV diagram needed to find a particular collinear limit, as follows, 

AM = 0 MHV: 1 + , 2 + , 3 + ,..., n + -> P + 

l-,2 + ,3 + ,...,n + ^p- 

AM = 1 ^ NMHV : 1~,2 + ,3 + ,... ,n + ^ P + 

1“,2",3 n + -> P~ 

AM = 2 NNMHV : l~,2~,3 + ,... ,n + ^ P + 

l - ,2",3",...,n + -> P~ 

(3.4) 

and so on for all AM > 2 cases. 

The splitting functions are derived by examining the general form of MHV dia¬ 
grams, which consist of MHV vertices and scalar propagators. The general form of 
the n-particle collinear splitting functions is given by 


> 

S* 

0 

"A ^ 

(3.5) 

AM 7 ^ 0 : oc - t-, 

' ■ v—1 ^ \^n—v 7 

(3.6) 


such that (v — 1) T (n — v) — n — 1, where v is the number of vertices, and thus v — 1 
is the number of scalar propagators. From (3.6) it follows that for an MHV-diagram 
to contribute to AM ^ 0 collinear limits, it is required to contain anti-holomorphic 
spinor products [ij\ of collinear momenta. However, because on-shell MHV vertices 
are entirely holomorphic, within the MHV rules there are only two potential sources 
of the anti-holomorphic spinor products. One source is scalar propagators 1 /s t] = 
1/ (i j) [j i] which inter-connect MHV vertices. The second source is the off-shell 
continuation of the corresponding connected legs in the MHV vertices. Each off-shell 
continued leg of momentum P gives rise to a factor (iP) oc (i|P1 77 ] which gives rise 
to anti-holomorphic factors of [jr)\. When the reference spinors rja are kept general, 
and specifically, not set to be equal to one of the momenta in the collinear set, the 
//-dependence must cancel and the off-shell continuation cannot give rise to an overall 
factor of [ij]. 

Therefore, the only source of singular anti-holomorphic factors are MHV-diagrams 
that contain an internal propagator of momentum = Pi+i + ... + Pj which is a 
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sum of external momenta from the collinear set such that q 2 = Si + ij —> 0. Hence, we 
conclude that only a subset of MHV-diagrams contributes to multi-collinear limits 
of tree amplitudes. The subset is determined by requiring that all v — 1 internal 
propagators are on-shell in the multi-collinear limit. This is a powerful constraint on 
the types of the contributing diagrams and it simplifies the calculation 3 . 

We exploit the universal nature of the splitting function by choosing to start 
with an amplitude with (n + 3) external legs, i.e. setting N = n + 3 in Eq. (3.3). 
The helicities of the gluons are adjusted so that the remnant ‘hard’ four point MHV 
amplitude 7L 4 (P A , (n + 1) + , (n + 2) -A , (n + 3) _ ) is given by 


H 4 ((n + 1) , (n + 2) , (n + 3) , P) — 

_ (n + 3, X) 4 _ 

(. P , n + 1) (n + 1, n + 2) (n + 2, n + 3) (n + 3, P) 


(3.7) 


with X = P for A = — and X = n + 2 for A = +. 

To read off the collinear limits from the MHV rules, we use the limiting expres¬ 
sions for the spinor products: (aq), ( bq) and ( ba ). Here a is a particle from the 
collinear set, b is a particle which is not in the collinear set, and q is the sum of the 
collinear momenta from i + 1 to j. Hence, using 


j j 

{aq} = (al) M , (bq) = Y ( bl > N > 

l = z ~|~1 l = i ~|~1 

and the expressions for spinors from the collinear set, 


(3.8) 


10 = Vzi\P) , If] = y/zi\P\ , |o) = sfz~a\P) , |o] = -JzAP] . (3.9) 

we have, 


j 


(aq) - 

* [ p v\ Y ( al )V^i = [ p vWi )(i,j;a) 

(3.10) 


l=i~h 1 


(bq) - 

3 

- IPriUbP) V z ‘ 

(3.11) 


l=i +1 


(ba) - 

(bP) \fz~ a . 

(3.12) 


Here we introduced the definition 

j 

A ( i)(i,j;a)= Y ( a 0 ( 3 - 13 ) 

1 

3 Note that this selection rule would not apply to neither gauge-fixed MHV rules (where r/’s are 
fixed to be equal to kinematic variables from the collinear set), nor to the BCF rules which mix 
holomorphic MHV vertices with anti-holomorphic MHV ones. 
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Equations (3.10) and (3.11) contain a factor [Pi]] which, however, will always 
cancel in expressions for relevant splitting functions. As such we can read off the 
collinear limits of the amplitudes from the MHV-rules expressions by replacing terms 
on the left hand side of equations (3.10), (3.11) and (3.12) with the expressions on 
the right hand side of those equations, and further dropping the [P rj\ factors. 

Certain terms in the sums that arise in MHV rules need special attention. These 
are the boundary terms involving either (0 1) or (nn + 1), and for these we have, 


{nn +1 ) _ 

A {1) (i,n;n+1) ^ EL:+i V 
( 01 ) y/zT 

A (1) (0,j;0) 


We now present our results for the splitting functions. 


(3.14) 

(3.15) 


4. Results 

In this section we give the results for the multiple collinear limit of gluons. First we 
give the general results for an arbitrary number of gluons with AM < 2. Afterwards 
we give explicit results for up to four collinear gluons for all independent helicity 
combinations, together with some specific examples for five and six collinear gluons. 

4.1 General results 

In this section we present the general results for the cases where the number of gluons 
with negative helicity changes by at most AM = 2, and those related by parity where 
the number of gluons with positive helicity changes at most by the same amount. 
With the help of parity these general splitting amplitudes are sufficient to obtain the 
explicit expressions for all helicity combinations of up to six gluons. 

We will often use a more compact notation for the splitting amplitude. We denote 
the splitting amplitude for n collinear gluons, of which r have negative helicity, by: 

split(l + ,..., ,..., m 2 , • • • j n^ri ■ ■ ■ j n+ P±) — Split±•'(mi,..., m r ). (4.1) 

4.1.1 AM = 0 

This is the simplest case which is read directly off the single MHV vertex. The 
denominator of an Appoint MHV amplitude is factorised as follows (in the limit of 
collinear pi,...,p n ): 

(N , 1) (1, 2)... (n, ?7 + 1)... (A — 1, N) 

n— 1 

= (^n^ + D) X ( (A/ - , P) (P, 77 + 1)... (AT — 1, AT) ^ (4.2) 



3+ 


Figure 2: MHY diagrams contributing to Split^(mi). Negative helicity gluons are indi¬ 
cated by solid lines, while arbitrary numbers of positive helicity gluons emitted from each 
vertex are shown as dotted arcs. 


where the first factor contributes to the splitting function, and the second one is the 
denominator of the remaining hard MHV amplitude. Hence, the splitting function 

is 


split (1 + ,..., n + —> P + ) 


1 

yj z l z n Ull ^ + 1 ) 


(4.3) 


and so by parity 


split (1 ,..., n —» P ) 


(-l)"" 1 

\/ z i z n nil ^+ 1 ] 


(4.4) 


Similarly, 


split (1 + ,..., m 1 ,..., n + —» P ) = 


mp 


\J z \ z n UlJ (l-i l + 1 ) 


(4.5) 


and 


split (1 — ,..., mf ,..., n~ —»• P + ) = mi 


\j z \ z n nr= i' [h i +1] 


(4.6) 


4.1.2 AM = 1 

This is the next-to-MHV (NMHV) case, and in the collinear limit we need to take 
into account only a subset of MHV diagrams. In fact, there is only a single MHV 
diagram (or more precisely a single class of MHV diagrams) which can contribute 
to Split(mi). It is shown in Fig. 2. 4 In the limit where gluons 1,... ,n become 
collinear. The left vertex in Fig. 2 produces a ‘hard’ MHV amplitude while the right 
vertex generates the splitting function. We need to sum over i and j in Fig. 2 in such 
a way that only diagrams with a singular propagator are selected in the collinear 
limit. This puts a constraint j < n where n is the number of collinear gluons. The 

4 MHV diagrams where hard negative helicity gluons are emitted from more than one vertex do 
not give rise to on-shell propagators and do not contribute in the singular limit. 
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Figure 3: MHV diagrams contributing to Split (mi, m 2 ). 


resulting splitting function reads, 


Split+(mi) = 


\/ziz n nL , 1 (i, 1 + 1 ) 


mi —1 n 

£ £ 

2=0 j=m\ 




(4.7) 


where we define 

D(i,j,q) = 


Qi+iJ 


{i, i + 1 ) (j, j + 1 ) 


A ( i)(i,i;i)A (1) (i,j;i + l)A (1) (i, j; j)A (1) (i, j; j + 1) 

(4.8) 


Similarly, there are three (classes of) MHV-diagrams contributing to Split (mi, m 2 ). 
They are shown in Fig. 3 and lead to a splitting function which reads 


Split^(mi,m 2 ) = 


mi — lm 2 — 1 2 a /• ■ \4 

^ A w^ , - 7;m i) 


m/ (1,1+ 1) v ^0 jt^i ^ 
S' E ^A ( 1 ) (i,i;m 2 ) 4 

i=m\ j=m 2 
mp —1 n 


££ 


(mi m 2 ) 


J? 9i+l 7) ' 

7=0 j=m 2 V + J' l=i +1 


£ 




(4.9) 


The remaining splitting amplitudes of the form 

split(1 _ ,..., m+,..., m -2 , • • •, m+,..., n _ —> P ± ) (4.10) 

are obtained by parity transformation through the usual replacement (/, /e) <->• — [/, h]. 
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4.1.3 AM = 2 

The collinear limits with AM = 2 are derived from next-to-next-to-MHV (NNMHV) 
diagrams. There are four (classes of) MHV-diagrams contributing to Split^(mi, m 2 ) 
which are shown in Fig. 4. The corresponding splitting function is, 

M) 


Splits (mi,m 2 ) = 


\jz\z n mu 1 + ^+1) 

e' E ”e E A « ( *’ j;mi)4A « (A; ’ r;m2)4 

i=0 j=m ,2 k=m\ r=m ,2 
m i — l fc m 2 —1 n 

SEES 

4=0 j=m\ k=m\ r=m 2 
/c n mi — 1 m 2 — l 

ESSE 

i=0 j=m 2 fc=0 r=m\ 
k n m \—1 j 


VD(i,j,q i+lJ ;k,r, q k+ i, r ) 

A(i)(f,i;mi) 4 A(i)(fc,r;m 2 ) 4 
DD(f, . 7 ,^+ 1 ,j-; fc, r, g*+i,r) 

A(i) (i, j ; m 2 ) 4 A ( i) (fc, r; mi ) 4 
DD(i,j,g i+ ij; /c,r, g fc+ i,r) 


SEES 


(mim 2 ) A (2 ) (z, j; k, r )' 
DD(i, 7, o i+ i k, r, qu+i r ) 

i=0 j=m 2 k =0 r=m 2 v 

where Ap)(i, j; k ) is given in Eq. (3.13) and we introduce 


(4.11) 


&{ 2 )(i,j;k,r) = E E Mv" 




41=4+1 l!=/c+l 

The ‘effective propagator’ DD is defined by 


(4.12) 


DD(i,j,„;fc,r, 92 ) = x('‘,k,qi,q 2 )x(r,i,q 2 ,qi)x(j,k,qi,q 2 )D(i,j,qi)D(k,r,q 2 ) 

(4.13) 


11 













Figure 5 : MHV topologies contributing to Split(mi, m 2 , m3). The negative helicity 
gluons mi, m 2 and m3 are distributed in a cyclic way around each diagram. The remaining 
leg is the negative helicity gluon that remains after the collinear limit is taken. 

in terms of D defined previously in Eq. (4.8), and x given by 

f 1 i 7^ k 

X(i, k, qi, q 2 ) = < A (2) ( qi , g2 )<«,i+i) . _ , • (4.14) 

[ A( 1 )(gi;i+l)A ( 1 )(g 2 ;i) 1 

Finally there are 16 classes of MHV-diagrams contributing to Split ^ (m 1 , m 2 , m 3 ), 
coming from the 5 topologies shown in Fig. 5 and their cyclic permutations. The 
individual contributions are given by 


Split^(mi, m2, m3) = 


16 


V^ULi (*> l + l )jz 


y^AW(mi,m 2 ,m 3 ) (4.15) 


where 


4 (1) (mi,m 2 ,m 3 ) 


A ( 2 ) (mi,m 2 ,m 3 ) 


m2 —1 n m3 — 1 j 

E E E E 4 

i=mi j=m,3 k=m2 r=ni3 


A(i) (i, j\ m 2 ) 4 A (1) (fc, r; m 3 ) 4 
DD(i, j,q i+ i,j-, k, r , q k +i, r ) 


k n m2 —1 m3 —1 

E E E E 4 

i=mi j=m 3 k=m\ r=m 2 


A(i) (i, j; m 3 ) 4 A (1) (fc, r; m 2 ) 4 
DD(bi,<h:+ij; A:,r, qfc+i.r) 
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mi — 1 m3 — 1 m2 — 1 j 


i (3) (mi,m 2 ,m 3 ) = EE EE EE EE 


2=0 j=rri2 k=m\ r=m2 
k m3 — 1 mi — 1 m2 — 1 


A { 4 \m ll m 2 ,rn 3 ) = EE EE 


1=0 j=m 2 k =0 r=m\ 
k n m 2 — 1 j 


i (5) (mi,m 2 ,m 3 ) = EE EE EE EE 


i=mi j=m 3 fc=mi r=ms 
k m3 — 1 mi — 1 j 


A (6) (mi,m 2 ,m 3 ) = JE E E 


A(i) (i, j; mi) 4 A ( i) (fc, r; m 2 ) 4 
13 DD(i, j,q i+ ij;k,r,q k+1>r ) 

A(i) (i, j ; m 2 ) 4 A ( i) (fc, r; mi ) 4 
1 DD(i,j,g i+ ij; A;, r, g fc+ i,r) 

(m 2 m 3 ) 4 A ( 2 ) (i, j; Ayr ) 4 
x DD(i, j,q i+ ij;k,r, q k + i, r ) 

(mim 2 ) A(2)(i,j) k,r ) 4 


A—- ms DD(i, j, g i+ i k, r, q k +i r ) 

i =0 j=m 2 k= 0 r=m 2 v 

mi -1 n m 3 -l j j , , \4 * / 7„ 


2 4 (7) (mi,m 2 ,m 3 ) = EEEE E 


i=0 j=m 3 k=m 2 r=ms l=i+ 1 
mi —1 n m2 —1 i j 


#(mi,m 2 ,m 3 ) = EE EE EE EE ( EE 


i =0 j=m 3 k=m\ r=m 3 Z=i +1 
k n mi — lm 3 — 1 j 


A ( - 9 \m 1 ,m 2 ,m 3 ) = EE EE EE EE ( EE 


2=0 j=ms k =0 r=m2 l=i +1 
mi —1 n m2 —1 m3 —1 j 


i (io) (mi,m 2 ,m 3 ) = EE EE EE EE ( EE 


i =0 j=m 3 k=mi r=m 2 Z=i +1 
Zc n mi — lm 2 —1 j 


A (11) (mi,m 2 ,m 3 ) = EEEE E 


2=0 j=ms k =0 r=m\ l=i+l 
m2 —1 n r m3 —1 a 


\ 4 (mi m 2 ) A ( i)(fc, r; m 3 ) 4 
/ DD(i,j,g i+ ij;/i;,r,gfc+i,r) 

\ 4 (m 2 m 3 ) 4 A ( i)(fc,r;mi) 4 
/ DD(i,j, k, r , (fc+i.r) 

4 (miffl 2 ) 4 A(i)(fc,r;m 3 ) 4 
DD(i, Ayr, g fc+ i,r) 

\ 4 K-m 3 ) 4 A ( i)(fc,r;m 2 ) 4 
/ DD(i, j, g i+ ij; k, r , q k +i, r ) 

4 (m 2 m 3 ) 4 A(i)(fc,r;mi) 4 
DD (i,j,q i+ 1 J ;k,r, q k+ i, r ) 


i (12) (mi,m 2 ,m 3 ) = EE EE EE EE 


i=m\ j=ms k=m2 r=m2 
mi — 1 m3 —1 r m2 —1 


J 4 (13) (mi,m 2 ,m 3 ) = EE EE EE EE 


A(1) (i, k ; m 2 ) 4 A(i) (r, j; m 3 ) 4 
VD(i,k,q i+lik ;r,j,q r+ 1 j) 

A ( i) (i, A;; mi) 4 A ( i) (r, j; m 2 ) z 


A^ A^ 4 —' A^ m s DD(i, A;, gj + i fcj r, j, q r+ \ b 

i= 0 j=m 2 k=mi r=m! v ’ ’ ) J?'ir+ijy 

mi —1 n r m 3 — 1 k , / \ 4 a ( 


A (14) (mym 2 ,m 3 ) = EEEE E 


2=0 j=ms k=m2 r=m2 l=i+ 1 
mi —1 n r m2 — 1 J 


# 5) (m 1 ,m 2 ,m 3 ) = EE EE EE EE ( EE 


2=0 j=ms fc=mi r=m\ l=r +1 
mi —1 n m2 —1 m3 — 1 a 


4 (m 1 m 2 ) A ( i)(r, j;m 3 ) 4 
DD(i, k, qi+i,k', r, j, q r +i,j) 

4 (m 2 m 3 ) 4 A(i)Q, fc;mi) 4 
DD (i,k,q i+ 1 y,r,j,q r+ltj ) 


A { 16 \m 1 ,m 2 ,m 3 ) = EE EE EE EE 


i =0 j=m,3 k=mi r=m 2 


A(i)Q, fc; mi) 4 A(i)(r, j\ m 3 ) 4 
VD(i,k,q i+li k;r,j,q r+ 1 j) 


4.2 Specific results for n < 7. 

In this section we present compact expressions for splitting amplitudes with up to 
six collinear gluons. These results are obtained directly from the general expressions 
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given in Section 4.1. 

First we note that splitting amplitudes satisfy reflection symmetry, 

split(l Al ,... ,n A " -»• P ± ) = (—l) n+1 split(n An ,..., l Al -> P ± ) (4.16) 

and the dual Ward identity, see e.g. [64], 

split(l Al , 2 A2 ,..., n Xn -»• P ± ) + split(2 A2 , l Al ,..., n A " -> P ± ) + • • • 

+ split(2 Az ,..., l Al , n An -»• P*) + split(2 A2 ,... ,n A ", l Al -> P ± ) = 0. (4.17) 

These relations reduce the number of independent splitting amplitudes significantly. 

4.2.1 n = 2 

For two collinear gluons there are two independent splitting amplitudes with AM = 
0. All others can be obtained by parity and reflection. Setting z\—z and z 2 = (1— z), 
we find 


split (1 + , 2+ - 

+ P + ) = 

1 

(4.18) 

a/ z( 1 — z) (12) 

split (1 — , 2 + - 

-P”) = 


(4.19) 

V4 1 - *) <12> 


As expected, the splitting amplitudes have a single pole proportional to (12) 1 . 
Note that in the soft limit z —> 0, we see that helicity conservation ensures that 
split (I - , 2+ -> P~) -> 0. 


4.2.2 n — 3 result from MHV rules 

For three collinear gluons there are three independent splitting amplitudes with 
AM = 0. They all follow directly from a single MHV vertex and are given by 


split (1 + , 2 + , 3 + - 

+ P + ) 

1 

y/Z\Z% (1 2) (2 3) 

split (1 — , 2 + , 3 + - 

+ P~) 

M 2 

■\Jz\Zz (1 2) (2 3) 

split(1 + , 2~, 3 + - 

+ P~) 

*2 

y/Z\Z$ (1 2) (2 3) 


Parity and the reflection symmetry, split (1 + , 2 + , 3 —» P ) = split (3 ,2 + ,l + 

P”), give the rest. 
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When AM = 1, there are three amplitudes, 


split (1 — , 2 + , 3 + -> P + ) = 


+ 


_ (1 2 ) ,- 2 2 _ 

*Jz\ZiZzS\2 (zi + z-i) ((1 3) ^Jz\ + (2 3) \fz-i) 
((12)y^+(13) \fz-zf 
Si, 3 (1 2) (2 3) ((1 3) v^T + (2 3) ’ 


(4.20) 


split(1 + , 2~, 3 + -*• P + ) = 


+ 

+ 


—split (2~, 1 + , 3 + -> P + ) - split (1 + , 3 + , 2~ -*• P + ) 

_ (12)^ 2 _ 

*Jz\ZiZzS \2 {z\ + z-i) ((1 3) yfz\ + (2 3) \fz-i) 

_ ((2 1) \fzj + (2 3) y/zj ) 4 _ 

3 1>3 (1 2) (2 3) ((1 3) y/zi + (2 3) ((1 2) ^ + (1 3) 

_MiL_ (4.21) 

\Z'^1^2^3'®2,3 (^2 + ^ 3 ) ((1 2) a/^2 + (1 3) \/53) 


split (1 + , 2 + , 3" -> P + ) = split (3“, 2 + , 1+ -»• P + ) . 


(4.22) 


In addition to singular terms like (12), we see that the splitting functions contain 
mixed terms like 3 ^ 3 . The net singularity is schematically of the form [ ]( ). 

Note that split(l _ , 2 + , 3 + —> P + ) contains poles in .s' 1.2 and the triple invariant 
S 13 = S 123 but not in s 2 , 3 - This is because there is no MHV rule graph with a 
three-point vertex involving two positive helicity gluons. 

Expressions for these splitting functions are given in Eq. (5.52) of Ref. [64], 
The results given here are more compact and have a rather different analytic form. 
After adjusting the normalisation of the colour matrices, the splitting functions of 
Eqs. (4.20)-(4.21) numerically agree with those of Ref. [64], 

4.2.3 n — 3 result from the BCF recursion relation 

We now want to rederive the above results using the BCF recursion relation of [3]. 
In doing this we will (a) draw some useful comparisons between the ‘BCF recursion’ 
and the ‘MHV rules’ formalisms from the perspective of collinear amplitudes; and 
(b) test our expressions, such as Eq. (4.20) for split(l _ , 2 + , 3 + —»• P + ). 

We start with the six-point amplitude A( l - , 2 + , 3 + , 4 + , 5“, 6 ”), and calculate 
it via the BCF recursive approach. We ultimately want to take the collinear limit 
1 j | 2 11 3 —> P + , so it will be convenient to choose the ‘marked’ gluons (required for the 
BCF recursive set-up) to be from this collinear set. Hence, we will mark the 1“ and 
2 + gluons. There are only two BCF diagrams which contribute to the full amplitude, 
and they are shown in Fig. 6 . We now note that in this particular collinear limit, only 
the second of these diagrams contains an on-shell propagator, l/s 2 3 - Nevertheless, 
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Figure 6: BCF diagrams contributing to A(1 ,2 + ,3 + ,4 + ,5 ,6 ). 


in distinction with the MHV rules approach which we have adopted previously, both 
BCF diagrams need to be taken into account in the collinear limit. 

The full amplitude reads 


A ( l “, 2 + , 3 + , 4 + , 5 _ , 6 ") = —-——— (-— + yt -1 


+ 


(1|2 + 3|4] ; 


( 3|1 + 216 ] \ [6 1][1 2](3 4)(4 5 ) s 3)5 [4 5 ] [5 6](1 2)(2 3 ) si , 3 

( 4 . 23 ) 

where the two terms on the right hand side correspond to the two BCF diagrams 
above (cf. Eq. (2.9) of Ref. [3]). 

In the 1 || 2 || 3 —» P + collinear limit, the first term becomes 


(1 2)4 


x 


(5 6) 4 


\/^i^2^3 s i,2(^i + ^2)((1 + (2 3 )^ 2 ) (5 6) (6 P)(P 4) (4 5) 


(4.24) 


This term factors into a contribution to the splitting amplitudes multiplied by a four- 
point MHV vertex. In contrast, in the collinear limit the second term factors onto 
the MHV type diagram, written in terms of the anti-holomorphic spinor products, 

((12)v^+(13)v^) 3 w [P 4] 4 


x 


S ,, 3 { 12 ) (2 3 )({ 13 ) ^ + (2 3 ) y/zfi [P 4 ] [4 5 ] [5 6][6 P]' 


( 4 . 25 ) 


For the special case of four-point amplitudes, the MHV and MHV amplitudes coincide 
and we find an identical result to Eq. (4.20). 

Likewise, to test our expression for split(1 + , 2~, 3 + —> P + ) we start from Eq. (3.4) 
in [3]; 


H ( l + , 2 -, 3 + , 4 -, 5 + , 6 -) = 


[1 3] 4 (4 6) 4 


[12] [2 3] <4 5 > <5 6> si j3 <6 11 + 213] <412 + 311] 
(2 6) 4 [3 5] 4 

(6 1 ) (1 2 ) [3 4 ] [4 5]s 3j5 (6|4 + 5 | 3 ] ( 2|3 + 4 | 5 ] 
[1 5] 4 (2 4) 4 

(2 3)(3 4)[5 6 ] [6 l ] s 2 ) 4 ( 4|2 + 3 | 1 ]( 2|3 + 4 | 5 ] 


. ( 4 . 26 ) 
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Taking the collinear limit 11| 2 || 3 —> P + , we find that 


Split (1 + , 2 ,3 


P + ) = 

+ 

+ 


441 1 2 ] 


\Jz\Z 2 Z 3 S 1 , 2(^1 + Z 2 ) ([1 3]yTi + [23]yT2) 

_[13T_ 

Si ,3 [1 2] [2 3] ([1 3] vTT + [2 3] ^ 2 ) ([1 2] ^ + [1 3] 

441 2 3] 


\Z^1^2^3 S 2,3(^2 + £ 3 ) ([1 2]y / i2 + [1 3]yT3) 


(4.27) 


This result has the same kinematic-invariant pole structure as Eq. (4.21), but other¬ 
wise is not obviously equivalent to Eq. (4.21). Note that Eq. (4.27) contains terms 
like ([1 2 ]y/z 2 + [1 3]yT3) (rather than ((1 c l) y /z 2 + (1 3)yT3)). Despite appearances, 
a more careful (e.g. numerical) comparison shows that these two results, Eqs. (4.21) 
and (4.27), are in fact the same. 


4.2.4 n — 4 

For n = 4, there are five collinear limits coming directly from MHV amplitudes where 
the number of gluons with negative hclicity doesn’t change, AM = 0, 


split (1 + , 2 + , 3 + 

, 4 + - 

-> P + ) 

1 

y/ziz 4 ( 12 ) (2 3) (3 4)’ 

(4.28) 

split (1 — , 2 + , 3 + 

,4 + - 

+ P~) 

4 

y/z\Z± (1 2) (2 3) (3 4)’ 

(4.29) 

split (1 + , 2 _ , 3 + 

,4 + - 

+ P~) 

4 

y/zizt{ 12) (2 3) (3 4) 

(4.30) 


The remaining two are obtained by reflection symmetry, 

split(1 + , 2 + , 3“, 4+ -> P-) = —split(4 + , 3", 2 + , 1+ -> P~), (4.31) 

split(1 + , 2+, 3+, 4- -> P~) = -split(4", 3 + , 2+, 1+ -> P~). (4.32) 

When AM = 1, there are ten splitting amplitudes however only three are imde- 
pendent, 


split(1 2+, 3 + , 4+ P + ) = £fi( 1, 2, 3,4) 

= _ 4' 2 ( 12 ) _ 

y/z 1 Z 4 (3 4) Si 2 {z\ + Z 2 ) A(i)(0, 2; 3) 

_ A(i)(0 ; 3;l) 3 _ 

\[zi (1 2) (2 3) s lj3 (zi + Z 2 + Z 3 ) A (1) (0, 3; 3)A(!)(0, 3; 4) 

A(!)(0,4; l) 3 

(12) (2 3) (3 4) Sl , 4 A (1) (0,4;4)’ 
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split(l - , 2~, 3 + , 4+ - P~) = B 2 { 1, 2, 3,4) 

_ z\ /2 zl /2 (2 3) _ 

V / ^2^4'S2,3A(i)(l, 3; 1)A(!)(1, 3; 4) 

_ ^i 3 / 2 A ( i)(l,4;2) 3 _ 

(2 3) (34) s 2 , 4 A (1) (1, 4; 1)A {1) (1, 4; 4) (1 - Zl ) 

_ (1 2 ) (z\ + z 2 ) 3 _ 

\JZiz 2 z 4 (34) Si,2A(i)(0, 2; 3) 

,_ (1 2) 3 (1 — Zj) 2, _ 

yfzl (2 3) S 1)3 A (1) (0, 3; 1)A (1) (0, 3; 3)A (1) (0, 3; 4) 


s M A (1) (0,4;l)A (1) (0,4;4) (23) (34)’ 


and, 


split(1", 2 + , 3",4+ ^ P~) = B 3 { 1, 2, 3,4) 

_ 4 7 y (i2) _ 

\fz 1^4 (3 4) Si )2 (zi + z 2 ) A {1) (0, 2; 3) 
zl /2 zl /2 (2 3) 

yj -2'3-2'4^2,3A(i)(1, 3; 1)A(!)(1, 3; 4) 

_ ^i 3/2 A ( i)(l, 4; 3) 4 _ 

(2 3) (3 4) S2, 4 A(i) (1,4; 1) A (1) (1,4; 2) A (1) (1, 4; 4) (1 - Zl ) 

_ ^i 3/2 A (1) (2,4;3) 3 _ 

(1 2) (3 4) s 3 , 4 A(i) (2,4; 2) A (1) (2,4; 4) (z 3 + z 4 ) 

(1 3) 4 (z\ + z 2 + Z 3 ) 3 

(1 2) (2 3) si j3 A (1 )( 0 , 3; 1)A (1) (0, 3; 3)A (1) (0, 3; 4) 


(1 2) (2 3) (3 4) s 1i4 A (1) (0, 4; 4) A (1) (0, 4; 1) ’ ' 

where Ais given in Eq. (3.13). The seven remaining AM = 1 splitting 
functions can be obtained by using the dual ward identity, 

split (1+, 2-, 3 + , 4+ P+) = -B\ (2,1, 3,4) - B 1 ( 2, 3,1,4) - B x { 2, 3,4,1), 

split(l + , 2 + , 3", 4+ -> P+) = B,( 3,4, 2, 1) + B 4 ( 3, 2, 4,1) + Bi(3, 2, 1,4), 

split (1~, 2+, 3 + , 4- -> P“) = B 3 ( 4, 3,1, 2) + B 2 (4, 1, 3, 2) + B 2 { 1,4, 3, 2), 

split(1 + , 2”, 3", 4+ - P~) = -B 3 ( 2,1, 3,4) - B 2 (2, 3,1,4) - B 2 ( 2, 3,4,1), 

(4.36) 


18 



or reflection symmetry, 


split(1 + , 2 + , 3 + , M - 

-»• P + ) 

= —split (4 — , 3 + , 2 + , 1+ - 

+ P~), 


split (1 + , 2“, 3 + , 4' - 

-+P-) 

= —split (4 — , 3 + , 2 _ , 1 + - 

+ P~), 


split(l + , 2 + , 3 _ , 4' - 

-+P-) 

= —split (4 — , 3 _ , 2 + , 1 + - 

+ P~)- 

(4.37) 


Finally splitting functions with AM = 2, 3 are related to those given above by 
the parity transformation. 

Inspection of Eqs. (4.33), (4.34) and (4.35) reveals that each term is inversely 
proportional to a single invariant, in keeping with its MHV rules origins. For this 
type of collinear limit, there are potentially six invariants, the double invariants 
S 12 , $ 2 , 3 , s 3 , 4 > the triple invariants si^, S 2,4 and Some poles are absent because the 
MHV rules forbid that type of contribution. For example, in split (1 — , 2 + , 3 + , 4 + —> 
P + ), there are no contributions with poles in s 2)3 , S 3,4 or s 2 ,4 precisely because these 
poles correspond to forbidden MHV diagrams. 

Expressions for the four gluon splitting functions are given in Ref. [64], The 
results given here are more compact and have a rather different analytic form. Af¬ 
ter adjusting the normalisation of the colour matrices, the splitting functions of 
Eqs. (4.33)-(4.35) numerically agree with those of Ref. [64], 

4.2.5 n = 5 

In total there are 64 different splitting amplitudes, but only eleven are independent. 
The rest can be obtained with the help of parity, reflection and dual ward identities. 
The three simplest independent collinear limits can be obtained using only MHV 
rules, 


split (1 + , 2 + , 3 + , 4 + , 5 + - 

■+ P + ) 

1 

y/ziz 5 ( 12) (2 3) (3 4) (4 5)’ 

split (1 — , 2 + , 3 + , 4 + , 5 + - 

+ P~) 

y/ziz 5 ( 12) (2 3) (3 4) (4 5)’ 

split(1 + , 2“, 3 + , 4 + , 5 + - 

+ P~) 

4 

y/ziz 5 ( 12) (2 3) (3 4) (4 5) 

The amplitudes with AM = 1 require the application of Eqs. (4.7) and (4.9). There 


are 5 independent amplitudes in this class of splitting amplitudes, but we give here 
only two examples, one for each of the cases- > + and-> —, 

split(1 — , 2 + , 3 + , 4 + , 5 + -> P + ) = 

__ (A ( 1 )(0,2;l)) 3 vTr _ 

v ^5 <1 2) (3 4) (4 5) s lj2 (Z! + z 2 ) A ( 1 ) (0, 2 ; 2)A ( 1 ) (0, 2 ; 3) 

_ (A ( 1 ) ( 0 , 3; l )) 3 yfz[ _ 

\fz 1 Z 5 (1 2 ) (2 3) (4 5) S 13 (zi + z 2 + Z 3 ) A ( 1 ) (0, 3; 3)A(!)(0, 3; 4) 
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(4.38) 


_ (A (1 )(0,4;l)) 3 y5r _ 

y/ziZ 5 (12) (2 3) (3 4) Si ; 4 (z\ + Z 2 + Z 3 + Z 4 ) A(i) (0,4; 4) A(i) (0,4; 5) 

_ (A,i)(0.5;l)) 3 

(12) (23) (34) (45)s,, 5 A (1) (0,5;5)’ 

and, 

split (1“, 2~, 3 + , 4 + , 5 + -*• P~) = 

_ z \ (A( 1) (l, 3; 2)) 3 _ 

(2 3) (4 5) s 2i3 A (1) (1, 3; 1) A (1) (1, 3; 3) A (1) (1, 3; 4) 

_ ^i 2 (A( 1) (l, 4; 2)) 3 _ 

(2 3) (3 4) s 2j4 A (1) (1, 4; 1) A (1) (1,4; 4) A (1) (1,4; 5) 

_ zi 2 (A (1) (l,5;2)) 3 ^ _ 

\fz\zl (2 3) (3 4) (4 5) s 2 ,5A(i)(l, 5; l)A(x)(l, 5; 5) (z 2 + Z 3 + Z 4 + ^ 5 ) 

_((12 )f( Zl + z 2 f^l _ 

(3 4) (4 5) si )2 A (1) (0, 2; 1) A (1) (0, 2; 2) A (1) (0, 2; 3) 

((1 2)) 3 (. z\ + Z 2 + ^ 3 ) 3 \fz{ 

(2 3) ( 4 5 > MsAq) ( 0 , 3; 1) A (1) (0, 3; 3) A (1) (0, 3; 4) 

((1 2)) 3 ( Z\ + Z 2 + Z 3 + ^ 4 ) 3 y/Zl 

7^ (2 3) (3 4) si, 4 A (1) (0, 4; 1)A (1) (0,4; 4)A (1) (0,4; 5) 


(2 3) (3 4) (4 5) s 1 )5 A (1) (0, 5; 1 )A ( 1 ) ( 0 , 5; 5) ’ ' 

The most complicated amplitudes are those with AM = 2 and require the use 
of Eq. (4.11). There are three independent splitting functions, but here we only give 
one example. 

split(1 2”, 3 + , 4 + , 5 + -*• P + ) = 

_ (A w (0,3; 1)) 3 (A w (1,3;2)) 3 _ 

\fzi (2 3) (4 5) A( 2 )( 0 , 3; 1 , 3 )si ; 3 (z\ + z-i + £ 3 ) A(p( 0 , 3; 4 )s 2 , 3 A(i)(l, 3; 1 )A(!)( 1 , 3; 3) 

_ (A ( 1 ) (0,4; l )) 3 (A( 1 } ( 1 , 3; 2)) 3 _ 

y/zE (2 3) s M (1 - z 5 ) A ( 1 ) ( 0 ,4; 4)A ( 1 ) (0,4; 5)s 2 ) 3 A ( 1 ) (l, 3; 1 )A ( 1 ) ( 1 , 3; 3)A ( 1 ) (1, 3; 4) 

_ (A ( 1 ) (0,4; 1 )) 3 (A ( 1 ) (l, 4; 2)) 3 _ 

(2 3) (3 4) A( 2 ) ( 0 ,4; 1 ,4)s M (1 - z 5 ) A (1) ( 0 ,4; 5)s 2 ) 4 A (1 ) ( 1 ,4; 1 ) A ( 1 ) (l, 4; 4) 

_ (A ( 1 ) (0, 5; l )) 3 (A ( 1 ) (l, 3; 2)) 3 _ 

A ( 1 ) (l, 3; 4) (2 3) (4 5) s ll 5 A ( 1 ) ( 0 , 5; 5)A ( 1 ) (1, 3; 1 )A ( 1 ) ( 1 , 3; 3)s 2 , 3 

_ (A ( 1 ) (l, 4; 2)) 3 (A ( 1 ) ( 0 ,5; l )) 3 _ 

s 2 ) 4 A (1) (1,4; 1 ) A (1) ( 1 ,4; 4) A (1) ( 1 ,4; 5) (2 3) (3 4) s 1 > 5 A (1) ( 0 , 5; 5) 
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(A ( 1 ,(1,5;2)) 3 (A ( 1 ) (0,5;1)) ! 


+ 


+ 


+ 


A( 2 ) ( 0 , 5; 1 , 5) A (1) ( 1 , 5; 1) A (1) (1, 5; 5)s 2 , 5 (2 3) (3 4) (4 5) s 1)5 

_ (( 12)) 3 (A ( 2 ) ( 0 , 3; 0 , 2)) 3 _ 

\fz 5 (45) Si ,3 (zi + z 2 + z 3 ) A ( 1 ) ( 0 , 3; 3)Ap)(0, 3; 4 )si j 2 A( 1 )(0, 2 ; l)Ap)( 0 , 2 ; 2 )Ap)( 0 , 2 ; 3) 

_ ((1 2)) 3 (A( 2 )( 0 ,4; 0 , 2)) 3 _ 

y/z 5 (3 4) Si j4 (1 — Z 5 ) A ( 1 ) ( 0 ,4; 4)Ap)(0,4; 5 )si ) 2 A(!)(0, 2 ; l)Ap)( 0 , 2 ; 2 )Ap)( 0 , 2 ; 3) 

_ ((1 2 )) 3 (A(2)(0, 4; 0, 3 )) 3 _ 

\fz 5 (2 3) s 4j4 (1 — Z 5 ) A(i)(0,4; 4)Ap)(0,4; 5 )si ) 3 A(i)( 0 , 3; l)A(i)( 0 , 3; 3)Ap)(0, 3; 4) 

_ ((1 2 )) 3 (A(2)(0, 5; 0, 2 )) 3 _ 

(3 4) (4 5) s 1 i 5 A (1) ( 0 , 5; 5)s 1 ) 2 A (1} ( 0 , 2 ; 1)A ( 1 ) (0, 2 ; 2 )A ( 1 ) ( 0 , 2 ; 3) 

_ (A(2)(0,5; 0,3)) 3 ((12)) 3 _ 

A ( 1 ) ( 0 , 3; 4)s 1 ) 3 A ( 1 ) (0, 3; 1 ) (2 3) (4 5) s 1 ) 5 A ( 1 ) ( 0 , 5; 5)A ( 1 ) (0, 3; 3) 

(A ( 2 )( 0 , 5; 0 ,4 )) 3 ((1 2)) 3 


A ( 1 ) ( 0 ,4; 5)A ( 1 ) (0,4; 1 ) (2 3) (3 4) s 1 ) 5 A ( 1 ) (0, 5; 5)s M A ( 1 ) (0,4; 4)' 


(4.40) 


4.2.6 n — 6 

Finally, for six collinear gluons there are 2 7 = 128 different splitting amplitudes, 
which can be expressed by 23 independent ones. To find all independent amplitudes 
we have to use Eq. (4.15) for the Erst time. Due to the length of the results we give 
here only two examples obtained with the help of Eqs. (4.7) and (4.9), 

split (1 _ , 2+, 3 + , 4 + , 5 + , 6 + -»• P + ) = 

(A { 1 ) ( 0 , 2 ; l )) 3 yfz[ 

(1 2) (3 4) (4 5) (5 6 ) si , 2 (^ + z 2 ) A (1) ( 0 , 2 ; 2 ) A (1) ( 0 , 2 ; 3) 

(A ( 1 ) ( 0 , 3; l )) 3 y/z[ 


+ 

+ 

+ 


y/zl ZE (1 2) (2 3) (4 5) (5 6 ) si )3 ( Zl + z 2 + z 3 ) A (1) (0, 3; 3) A (1) (0, 3; 4) 

_ (A ( 1 ) (0,4;l)) 3 y7I _ 

y/z\Z 6 ( 12 ) (2 3) (3 4) (5 6 ) Si )4 (zi + z 2 + z 3 + Z 4 ) Ap)( 0 ,4; 4)A(i)(0,4; 5) 

_ (A ( 1 ) (0, 5; l )) 3 y/zi _ 

yfziz 6 ( 12 ) (2 3) (3 4) (4 5) 3^5 (z 3 + z 2 + z 3 + Z 4 + z$) Ap^O, 5; 5)Ap)(0, 5; 6 ) 
(A ( 1 ) ( 0 , 6 ; l )) 3 


(1 2 ) (2 3) (3 4) (4 5) (5 6 ) s li 6 A (1) ( 0 , 6 ; 6 ) ’ 


(4.41) 


split(1 _ , 2", 3 + , 4 + , 5 + , 6 + -»• P~) = 

Z\ (A ( 1 ) (l, 3; 2)) 3 


(2 3) (4 5) (5 6 ) (1, 3; 1 ) A ( 1 ) (l, 3; 3) A ( 1 ) ( 1 , 3; 4) 
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_ ^i 2 (A (1 )(l,4;2)) 3 _ 

\/ziZ6 (2 3) (3 4) (5 6) s 2 , 4 A (1) (1,4; 1) A (1) (l, 4; 4) A (1) (1,4; 5) 

_ z i 2 (A(!)(l, 5; 2)) 3 _ 

(2 3) (3 4) (4 5) s 2 , 5 A (1) (1, 5; 1)A (1) (l, 5; 5)A (1) (1, 5; 6) 

_ ^i 2 (A (1) (l,6;2)) 3 y^ _ 

\/z\Z§ (2 3) (3 4) (4 5) (5 6) s 2 ,6A(i)(l, 6; l)A(i)(l, 6; 6) (z 2 + z^ + Z4 + z$ + Zq) 

_ ((i2)) 3 (^ + ^ 2 ) 3 ^r _ 

(3 4) (4 5) (5 6) 5 1>2 A (1) (0, 2; 1) A (1) (0, 2; 2) A (1) (0, 2; 3) 

((1 2)) 3 (zi + z 2 + ^3) 3 y/zi 
(2 3) (4 5) (5 6) s 1i3 A (1) (0, 3; 1) A (1) (0, 3; 3) A (1) (0, 3; 4) 

((1 2)) 3 (zi + Z2 + Z3 + 24) 3 *Jz[ 

^z % (2 3) (3 4) (5 6) s li4 A (1) (0,4; 1) A (1) (0,4; 4) A (1) (0,4; 5) 


((1 2)) 3 (zi + Z2 + Z3 + Z4 + ^ 5 ) 3 yfz[ 

^z % (2 3) (3 4) (4 5) s 1>5 A (1) (0, 5; 1)A (1) (0, 5; 5)A (1) (0, 5; 6) 

_ (( 12)) 3 _ 

(2 3) (3 4) (4 5) (5 6) s 1)6 A (1) (0, 6; 1) A (1) (0, 6; 6)' 


(4.42) 


5. Conclusion 

In this paper we have considered the collinear limit of multi-gluon QCD amplitudes 
at tree level. We have used the new MHV rules for constructing colour ordered am¬ 
plitudes from MHV vertices together with the general collinear factorization formula 
to derive timclike splitting functions that are valid for specific numbers of negative 
hclicity gluons with an arbitrary number of positive helicity gluons (or vice versa). 
In this limit, the full amplitude factorises into an MHV vertex multiplied by a multi- 
collinear splitting function that depends on the helicities of the collinear gluons. 
These splitting functions are derived directly using MHV rules. Out of the full set of 
MHV-diagrams contributing to the full amplitude, only the subset of MHV-diagrams 
which contain an internal propagator which goes on-shell in the multi-collincar limit 
contribute. 

We find that the splitting functions can be characterised by AM, the difference 
between the number of negative helicity gluons before taking the collinear limit, and 
the number after. AM + 1 also coincides with the number of MHV vertices involved 
in the splitting functions. Our main results are splitting functions for arbitrary 
numbers of gluons where AM = 0,1,2. Splitting functions where the difference 
in the number of positive hclicity gluons A P = 0,1,2 are obtained by the parity 
transformation. These general results are sufficient to describe all collinear limits 
with up to six gluons. We have given explicit results for up to four collinear gluons 
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for all independent hclicity combinations, which numerically agree with the results 
of Ref. [64], together with new results for five and six collinear gluons. This method 
could be applied to higher numbers of negative helicity gluons, and via the MHV- 
rules for quark vertices, to the collinear limits of quarks and gluons. 

We anticipate that the results presented here will be useful in developing higher 
order perturbative predictions for observable quantities, such as jet cross sections at 
the LHC. 
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